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Abstract. We discuss the phenomenology of the minimal inflation scenario. We concentrate on two 
aspects: inflationary trajectories and particle production. Our findings can be summarized in two 
main results: first, that inflationary trayectories that are very fiat and provide enough number of 
e-foldings are natural in the scenario without fine tunning. We present a general formalism to identify 
attractors in multi-field infiation regardless of trajectories fulfilling the slow-roll conditions. We then 
explore particle production in the model and show how the infiaton naturally transmutes into a dark 
matter particle. One interesting feature of our model is that it provides a novel mechanism to generate 
particles and entropy in the universe: the filling of the Fermi sphere up to a given momentum pp due 
to the sea of goldstinos that are an important part of the matter generated after infiation. With this 
mechanism in hand we predict that the gravitino should have a mass of > 100 — 1000 TeV. Another 
interesting feature of our model is that the predicted level of gravity waves is ?' = 0.1 — 0.001, which 
is in the range of detectability from Planck and upcoming CMB polarization experiments. 
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1 Introduction 

In recent realizations of the theory of inflation[l], it is most natural to use niulti-field inflation, that is 
the inflaton field is a composite of many fields, with some recent examples where the number of fields 
can be very large [2] . It is difficult to think of examples of the inflaton based on fundamental physics 
models where one can be limited to only a single field. While it can be the case that many of the fields 
are "frozen" and that only one of them drives the inflationary period, it is most often the case that 
several fields participate in the inflationary period. In order to determine the existence of attractors 
in multi-field inflation it is common practice, in order to determine if a given trajectory produces 
enough number of e-foldings to explain the size of the current universe, to impose the so-called slow- 
roll condition simultaneously in all fields that drive inflation [3]. This is, of course, a sufficient but not 
a necessary condition. It remains an open problem in multi-field inflation to determine the attractor 
in phase-space. Further, in multifield infiation the infiaton trajectories can be dominated in part by 
non-linear friction terms. Motivated by an interest in the existence of the attractor trajectories in 
multi-field inflation and to study the phenomenology of the minimal infiation scenario [4, 5] we have 
decided to revisit the definition of "slow-roll" conditions in multi-field models of infiation. 

We have also investigated particle production in the minimal infiation scenario in order to explore 
the fate of the inflaton. The reason we have done this, and not used standard (pre)(re)-heating calcu- 
lations for a scalar field, is that in some minimal infiation scenarios the inflaton is not a fundamental 
field but a bilinear fermion, or a combination of a fermion bilinear and a scalar field. Therefore, under 
certain conditions the infiaton field can "transmute" in part into a sea of fermions and end up as part 
of the constituents of the Universe. 

In the minimal-inflation[4, 5] scenario the field X that drives the exponential expansion of the 
Universe can often be represented at low energies by a Goldstino composite GG. Our main motivation 
to propose to identify the infiaton field with the order parameter of supersymmetry breaking is guided 
by the fact that, independently of the particular microscopic mechanism driving supersymmetry 
breaking (in what follows we will restrict ourselves to i^-breaking ) we can define a superfield X 
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whose component at large distances becomes the "Goldstino" (see [6, 7]). In the UV the scalar 
component x of X is well defined as a fundamental field while in the IR, once supersymmetry is 
spontaneously broken, this scalar field may be expressed as a two Goldstino state. The explicit 
realisation of a; as a fermion bilinear depends on the low-energy details of the model. In models of 
low-energy supersymmetry the realization of x as GG can be implemented by imposing a non linear 
constraint in the IR for the X field of the type = 0. In our previous approach to inflation we 
used one real component of the UV x field as the inflaton. We assumed the existence of a F-breaking 
effective superpotential for the X-superfield and we induced a potential for x from gravitational 
corrections to the Kahler potential. In the cosmological context we have to be careful to verify if the 
conditions for the decoupling of the scalar field x, or rather, its transmutation into a goldstino bilinear 
are realised after the system exits infiation and begins the pre-heating period. In the simplest models 
analysed below this condition is not completely satisfied, hence, after inflation, we obtain a universe 
flUed with a fermi gas (liquid) generated by the goldstinos together with a massive scalar component. 

It is this property of this model of inflaton that allows us to use the inflaton as the generator 
of dark matter. When the inflaton abandons the slow- roll conditions (77 > 1; where r] is the second 
slow-roll parameter), we can work out the pre-heating and reheating of the universe. We have thus 
a possible candidate for dark matter. Since our model only depends on the scale of supersymmetry 
breaking, we get tight constraints for the expected mass of the dark matter particles 100 TeV). 



In multi-fleld inflation the question of finding the attractor for inflationary trajectories is more com- 
plicated than with a single field. In the case of caothic inflation the attractor can be found analytically 
[8]. Here we describe a general algorithm to obtain attractors in multi-fleld inflation. 

With a flat FLRW ansatz for the space-time metric, we consider the equations of motion for a 
generic model. The scalar fields will take values on some target space manifold with metric (7ij($). 
The first FLRW equation relates the Hubble parameter to the energy density: 



where the are the scalar fields of the theory, with the inflaton among them, and we have the usual 
deflnitions for pressure (p) and energy density (p) 
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(2.1) 



p = \g,,¥^^ - y($); p = \g,j¥^^ + V{^). 



(2.2) 



We also have the second cosmological equation: 



a 



AttG 



(2.3) 



a 
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To have inflation, we require that p = — p, as this will give an exponential growth of the scale 
factor a = exp(/ H{t)dt). Thus we need to find trajectories with p ^ —p for enough time: more than 
55 e-folds to explain the current universe. The equation of motion for the scalar field reads: 



dt 



3H^' + g'^Vji^) = (2.4) 



where D is the covariant derivative with respect to the target space metric Qij. 

In our case, the scalar fields form a complex scalar field, the partner of the goldstino field. 
Although we will make explicit the arguments below in this case, the arguments easily generalise to 
a more complex situation. Our complex field can be written as 2; = M{a + z/3)/-\/2 (hereafter M is 
the Planck mass scale). The potential is V{z, z) = f^V{a, /3), since we only include for simplicity two 
scales, the Planck scale M and the supersymmetry breaking scale /^/^. In supergravity models, the 
gravitino mass is up to simple numerical factors given by 7713/2 ~ f /M. It is convenient to write down 
dimensionless equations of motion such that time is counted in units of f~^^^. If is the space-time 
volume, then the Kahler metric reads 

ds'^ = 2g^^dsdz = d^d^K{a, P)M'^{da'^ + dl3^). (2.5) 



The action is then: 



S = L'J dta^ ^^(a, /3)M2(d2 + _ f2y^^^ ^ (2.6) 



where t = t Mj f, thus 

S = L^fm-^, J dra' (^^g{a, ^){a'^ + P'^) - V{a, /?)) . (2.7) 

Hence the system of differential equations for the trayectory becomes: 

a" + 3-a' + loggia'' - (3'^) + dp logga'p' + = 

p" + 3-/3' + Idff logff(/3'2 - + log ga'p' + g-% = 

CI ^ 

" ^ f-a{a'' + n + V{a,P)) (2.8) 



a m3/2 \/3 V2"' 

Next we need to understand the conditions leading to slow roll, and also the corresponding basin of 
attraction in the space of initial conditions leading to them. In general there is no need to start with 
conditions where the scalar field is at rest. They will evolve in time to get closer and closer to the 
slow roll attractor. To determine the attractor we notice that the slow roll conditions imply that in 
(2.4) we can ignore the acceleration term D^^ /dt which represents geodesies in the target space. 
This implies that the trajectory is dominated by the lower derivatives terms a' = —g~^Va/{3H), 13' = 
—g~^Vp/{3H). The initial conditions for the trajectories are then provided by these equations at 
some fiducial initial time which yields q:'(0), /3'(0), but with the Hubble parameter H determined by 
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(2.1). Replacing the derivatives determined by the first order equations in (2.1) yields a fourth order 
equation for H. The only acceptable solution is: 



^ {iV + ^/&V' + 9^2^ (2.9) 



^ -'18 

We now have the relevant initial conditions for the geodesic equation D^^ /dt ^ 0. We start with the 
same initial conditions for the geodesic equation and the first order equations representing slow roll, 
and analyse how long the trajectories stay together. This is shown in the right panel of Fig. I. So 
long as the geodesies and the exact solution to the equations of motion remain close enough, the slow 
roll conditions hold. This argument should hold for the general case. Rather than looking for general 
analytic criteria, we now look for a numerical sampling of the initial conditions to find the attractor 
trajectories. To illustrate the existence of inflationary trajectories and attractors in phase space for 
multi-field inflation we focus on two examples in our minimal inflation scenario: the massles case and 
the "half-tube" case to be explained below. 

3 Some examples of minimal inflation potentials with very flat directions 

In our minimal inflationary scenario [4, 5] we use only the Ferrara-Zumino (FZ) multiplet to drive 
inflation. The scalar potential in the Eisntein frame is given by: 

o _ 

VE = e^{-j^WW + G^'^DxWDxW), (3.1) 

where the Kahler metric and the Kahler covariant derivatives are given by: 

Gxx ^ dxdxK{X,X) DW{X) ^ dxW{X) + ^ dx KW{X) (3.2) 

In our approach we make an explicit, but resaonably generic choice for K and W. 

For us the inflaton superfield is the FZ-chiral superfield X — z+y/2 9ijj+9'^F, the order parameter 
of supersymmetry breaking. We will consider the simplest superpotential implementing F-breaking 
of supersymmetry. More elaborate superpotentials often reduce to this one once heavy fields are 
integrated out. 

W = fX + foM (3.3) 

with /o some constant to be fixed later by imposing the existence of a global minimum with vanishing 
cosmological constant and with / the supersymmetry breaking scale / ~ iJ'lusy 

We are interested in subplanckian inflation, and not on the ultraviolet complete theory that 
should underlie the scenario. Hence we simply parametrize the subplanckian theory in terms of the 
previous superpotential, and a general Kahler potential whose coefficient will be taken of order one. 
We try to use our ignorance of the ultraviolet theory to our advantage. The Kahler potential K we 
consider is: 
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The most straightforward interpretation of why we write this Kahler potential is that we consider 
aU terms up-to power and a log term explicitly breaking all possible R-symmetries. The 

coefficients /o,a, fo, c will be chosen appropriately so that we obtain flat directions. For our purposes 
it is convenient to chose a = 0, as this guarantees the existence of a global minimum, fo will be 
adjusted such that the global minimum is at a vanishing value of the potential. So the model only 
has 6, c as free parameters, which we try to keep of order one to avoid fine tunning in the potential. 

From the collection of potentials considered, not all will show flat directions where it is possible to 
inflate during enough of e-foldings (> 55) for any choice of the microscopic parameters fo, a, b, c. We 
will restrict to cases where the potential has a global minimum with vanishing cosmological constant 
and thus we fix the value of the minimum at tuning the value of /q. 

For simplicity consider first the massless case, where /o = — /, a — 0,b ~ I, c — 0. This is 
the simplest minimal model as it has only one parameter 6, which we set to 1 for concreteness. The 
potential has the shape shown in the left panel of Fig. I. The obvious flat direction runs parallel to 
the (3 direction. If we solve the equations of motion, the slow-roll trajectories look like those shown 
in the right panel of Fig. I, where we have chosen arbitrary points to start the fleld. Note that all 
trajectories flow toward an attractor parallel to the /3 direction, i.e. the one that is quite flat and 
therefore the one providing enough e-foldings. 

A very interesting case of potentials with very flat directions are those where a — and c — ~g x 6, 
where g is a coefficient order one. This provides a very flat direction along a "half-pipe" with the 
added advantage that the steep directions provide a large set of initial conditions leading to inflation. 
In these models a global minimum always exits. The value of the potential at this minimum is not 
necessarily zero by simply choosing /o = — /, so fo needs to be chosen differently. To do this let us 
define fo — ^hf. Since the potential is symmetric in (3 and the global minimum, for these models, is 
always at /3 = 0, we only need to minimise the potential for values of h, g, a. Because g is order one in 
all cases, the potential will depend only on h, a, or rather the dependence on g is very weak. Therefore, 
the equation to find h,a will be V{h,a) — (see left-bottom panel of Fig. II). This equation has 
many solutions, but only the one for which fj^ = will provide a solution with a global minimum 
equal to zero. This can be done by solving the following implicit equation (where the dependence in 
g has been dropped because it is very weak) 



A fitting formula providing an accurate solution is /o = — 1 + 1-2 x lO^^g. This provides the full 
solution to find values for /o, 6, c (recall that we always keep a — 0) where the global minimum is at 
and a very flat direction exist. Let us study some examples. 




(3.5) 
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Figure I. Left-panel: the potential (shown in units of V/f^) for the massless case as a function of the real 
fields Q, p. Right-panel: inflationary trajectories computed with the formalism described in the text to choose 
initial conditions. Note the existence of an attractor that brings all the trajectories to the flat region of the 
potential. 

If we choose a — 0,b — l,c — — 1.7, /o = —.99819156/ then we obtain the potential shown in 
the upper-left panel of Fig. II. This potential has a very flat direction for sub-planckian values in 
both directions a and /3. This can be easily checked by computing the first slow-roll parameter in 
both directions and adding them e = g^-'V, iV,j/V^ (upper-right panel of Fig. II). This is of course a 
sufficient but not necessary condition for infiation. Wherever this number is smaller than one there is 
inflation (see also Fig. III). 

We can now compute trajectories and attractor in more detail using the equations from the 
previous section. This is shown in the lower-right panel of Fig. II. As expected all trajectories pass 
through the flat region, therefore any inflationary trajectory provides sufficient number of e-foldings. 
The bottom left panel of Fig. II shows the value of the potential along the inflationary trajectory as 
a function of time. 
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Figure II. Upper-left panel: The potential (shown in units of V/f^) as a function of the real fields (a, /?) for 
the case with parameters (a = 0,b = 1, c — —1.7, fo ~ — /, note that this creates a very flat " bottom- valley" 
region, where slow-roll conditions are fulfilled. The fact that this region is very flat is shown in the upper-right 
panel, where we show the value of the first slow-roll parameter added in the two real directions. Trajectories 
for arbitrary initial conditions, as chosen by the algorithm described in the text, are shown in the middle-left 
panel. Note that all converge to an attractor around the flat region. Middle-right panel; the minimisation of 
/i as a function of a in order to determine the minimum of the potential with zero value (or to match it to 
the present value of the cosmological constant). Bottom-left panel: the value of the potential as a function 
of time t (as defined in text) for the blue trajectory; the flat region is clearly visible as well as the exit from 
slow-roll, which starts particle production, yielding a canonical inflationary behavior 
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Figure III. The traditional slow-roll parameters (e, rj) are computed for different values of the potential in 
the "half-pipe" potential show that in all cases both e, 77 < 1. The values of the potential are chosen randomly 
from a monte-carlo realization of 1000 trials that keeps those values that fuUfil the slow-roll conditions. The 
values, from left-to-right, are: {a — 0,b — l,c — —1.7), (a = 0, 6 = 1, c = —1.65), {a = 0,b = 1.1, c = —1.65), 
(a = 0,6 = 1.15, c= -1.7). 

The level of gravitational waves is easily predicted in our model. Recall that r « (V/3.3. x 
lO^^GeV)'*. In our case, for inflationary trajectories that reproduce the amplitude of the fluctuations 
of the CMB, V - p, where is in the range lO^-^^ GeV. Therefore r - 0.1 - 0.001, which is in 
the range of detectability from Planck and upcoming CMB polarization experiments. 

4 Particle production at the end of inflation. 

After inflation ends, energy and entropy in the form of radiation need to be produced in order to 
describe the observed universe. One interesting feature of our choice of the field X as the inflaton, is 
the fact that because it is a bilinear fermion, at low energies it will manifest itself as a fermion rather 
than a scalar. There is then the possibility that gravitino production can be achieved directly without 
the need of parametric reheating. This is the scenario we explore in this section. 

When supersymmetry is spontaneously broken, we can follow the flow of X to the infrared (IR). 
In the deep IR this field satisfies a nonlinear constraint and becomes the goldstino superfield. 




(4.1) 



G2 



+ V2eG + e'^F 



(4.2) 



2F 
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The scalar component x of X becomes a goldstino bilinear. Its fermionic component is the 
goldstino fermion G, and F is the auxiliary field that gets the vacuum expectation value. This 
nonlinear constraint is a consequence of integrating out the x-field. To do this, the effective inflaton 
mass needs to be larger than the mass of the gravitino. To be more precise, toinf > 2m3/2 would be 
the required condition for the light fields to be integrated out. When we talk about integrating out 
we always refer to the above non- linear constraint being realised (see e.g. [7]). Ref. [7] have shown 
that this non-linear constraint occurs always in Minkowski space. Although, after the inflationary 
period one is in curved space (FLRW) , it is not difficult to see that corrections to the Minkowski result 
{X'^ — 0) will come as powers of 1/AP, and thus corrections will be suppressed for sub-plankian fields, 
as is the case of our model. The masses of the inflaton and the gravitino depend on the parameters of 
the model. In minimal models like the ones considered in this paper there is no room to have complete 
decoupling of the scalar field, which is only realised partially, hence the final state after inflation will 
contain a mixture of a Fermi liquid state in fermions together with a gas of inflatons. 

Let us compute the masses of the inflaton in the real vacuum for the "half-pipe" model. The 
masses in the two real directions a, /3 can be well approximated by the following formula 

miNF = 0.5^ (4.3) 

This is very similar to the gravitino mass 7713/2 — f/M, which is the mass obtained for the 
gravitino in most models of supersymmetry breaking. We can be more precise and compute the 
gravitino mass as 7713/2 = \W\e^^ /\/3. For several values of the coefficients a, b, c this gives values of 
m3/2 = (0.2 — 0.5)-^. So this reinforces the fact that both masses are very similar and that the non- 
linear constraint is not fully realised. Conversion of X into GG will take place as soon as the inflaton 
fleld abandons the slow-roll condition, where its effective mass is zero, and adquires an effective mass. 
If we assume that the most favorable scenario takes place we can immediately see that the Fermi 
sphere will be flUed up to momentum pf — \J -{j'fnwF- We can now estimate a lower bound for the 
gravitino mass to be able to produce enough of them at a rate that overcomes the expansion of the 
universe-recall that we arc now in a FLRW phase. The observed number of particles in the universe is 
^ 10™^^". Using the computed value for the Fermi momentum, this implies that y/f > lO^^GeV 
which gives a lower bound on the gravitino mass: m-^ > 100 — 1000 TeV. 

However, as we have seen above, because the similarities of the mass of both the inflaton and 
goldstino fields, it does not seem to be a very efficient route of gravitino production, even though some 
conversion will take place. We have then to explore in more conventional terms particle production. 
This is done following standard arguments, and we will not dwell on them here. We will come back 
to a more realistic treatment (including creation or ordinary matter, not just gravitinos and inflatons 
in a future publication). The basic mechanism to produce particles after inflation is to use broad 
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parametric resonances via oscillations of the fields around its minimum (see for instance Ref. [8] and 
references therein). For the fermions, we can follow the arguments in Ref. [9, 10]. 

The simplest way to proceed is to first choose the unitary gauge in the full supergravity la- 
grangian, and given that we are working below the Planck scale, we can safely ignore the quartic 
gravitino coupling. The only quartic terms that one should keep are those associated to the sigma 
model associated to the rigid supersymmetric sigma-model representing the naive M — > cxi limit. The 
relevant coupling of the X-field to the gravitino is through the term We^M^ if^a'^'^ipb/y/S. Although 
the problem seems complicated, it can be shown [9] that we can decompose the problem in two parts: 
the production of the spin-3/2 components of the gravitino that is highly suppressed, and the produc- 
tion of the spin- 1/2 component (effectively the absorbed goldstino given the gravitino a mass), which 
is not suppressed by powers of M. Furthermore, after some manipulations the authors of Ref. [9] show 
that the problem is formally identical with the production of ordinary spin-1/2 fermions coupled to 
the inflaton field. The lagrangian takes the form: 

^^^A72A/^ /o + /^ ^^^ (4.4) 

for all practical purposes the contribution of the K-prefactor is negligible. We can now use results 
from the literature to compute the particulc production rate [11-14]. In brief, they quantize the field 
following the standard semi-classical approximation to quantize fermions (Bogolyubov which in turn 
leads to a Hartree-Fock eq.). The important result is that the eigenmode solution obeys an oscillator- 
like equation with complex frequency, then one proceeds to solve this equation to find the particle 
density in an expanding background. In the case of a flat universe, we can separate variables, each 
mode is characterised by a given momentum k, and to a first approximation we end up with a Fermi 
sphere in FLRW. The interesting result in Ref. [14] is that only very few (2-5) oscillations produce 
values of of order a few tenths (remember that because we are dealing with fermions has to be 
below 1/2). Even for just two oscillations ^ 0.1. One can see that the Fermi sphere fills rather 
quickly, with pp ~ \J jjmi^-p. Which requiring again that the Fermi sphere fills to a level such that 
we produce enough particles (~ 10™~^°) gives a limit on the mass of the gravitino > 1000 — 10000 
TeV. Here we are not considering the contribution to the matter content of the inflaton field. In 
the type of models presented, the X-field does not fully decouple. The ideal scenario would be the 
one where at the end of inflation, the inflaton "integrates" itself out. this would require its mass 
to be sufficiently higher that the gravitino mass, so that it can be safely integrate out, and then we 
enter the non-linear regime where soon after inflation we have a Fermi gas of goldstinos in the FLRW 
universe. This of course requires considering more elaborate models, where one should also include 
the possibility of producing ordinary matter (squarks, quarks, leptons etc) that we have so far ignored 
systematically. 

There are however some results from the literature that hint toward the fact that baryons can be 
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produced in the model (see e.g. the review in Ref. [15]). To be more specific, there is however a couple 
of routes one could explore to achieve the desire production of baryons, or more precisely to excite 
solely the MSSM degrees of freedom: one could argue that inflation always ends in a vacuum where all 
the minimal super-symetric model (MSSM) gauge group is recovered; one advantage of doing so is to 
embed inflation completely within a visible sector, i.e. MSSM inflation models where inflation occurs 
via gauge invariant F-flat directions. Alternatively, one would have to know how inflaton couples to 
all the hidden sectors and the visible sector. In certain examples, such as in the case of large volume 
compactiflcations, where due to hierarchy in scales and the fact that SUSY is broken by fluxes, one 
can actually minimize the hidden sectors to some extent. We hope to come back to these issues in a 
future publication. 

Finally, we want to conclude remarking the fact that one firm prediction of our model is that 
the SUSY breaking scale has to be in the range 10^'^"^^ GeV, which in turns gives large masses for 
the gravitino (using that rrigiav = f /M). In order to determine the fate of the gravitino as a possible 
component of the dark matter we will investigate in a future publication how our model interacts with 
the MSSM (or similar realizations) to also explore how baryons are being produced. This requires 
knowledge of the microscopic details. 
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